Low mutation rate of spontaneous mutants enables detection of causative genes by comparing whole genome sequences.

In the early 1900s, mutation breeding to select varieties with desirable traits using spontaneous mutation was actively conducted around the world, including Japan. In rice, the number of fixed mutations per generation was estimated to be 1.38-2.25. Although this low mutation rate was a major problem for breeding in those days, in the modern era with the development of next-generation sequencing (NGS) technology, it was conversely considered to be an advantage for efficient gene identification. In this paper, we proposed an in silico approach using NGS to compare the whole genome sequence of a spontaneous mutant with that of a closely related strain with a nearly identical genome, to find polymorphisms that differ between them, and to identify the causal gene by predicting the functional variation of the gene caused by the polymorphism. Using this approach, we found four causal genes for the dwarf mutation, the round shape grain mutation and the awnless mutation. Three of these genes were the same as those previously reported, but one was a novel gene involved in awn formation. The novel gene was isolated from Bozu-Aikoku, a mutant of Aikoku with the awnless trait, in which nine polymorphisms were predicted to alter gene function by their whole-genome comparison. Based on the information on gene function and tissue-specific expression patterns of these candidate genes, Os03g0115700/LOC_Os03g02460, annotated as a short-chain dehydrogenase/reductase SDR family protein, is most likely to be involved in the awnless mutation. Indeed, complementation tests by transformation showed that it is involved in awn formation. Thus, this method is an effective way to accelerate genome breeding of various crop species by enabling the identification of useful genes that can be used for crop breeding with minimal effort for NGS analysis.

Genome-wide association study of leaf photosynthesis using a high-throughput gas exchange system in rice.

Enhancing leaf photosynthetic capacity is essential for improving the yield of rice (Oryza sativa L.). Although the exploitation of natural genetic resources is considered a promising approach to enhance photosynthetic capacity, genomic factors related to the genetic diversity of leaf photosynthetic capacity have yet to be fully elucidated due to the limitation of measurement efficiency. In this study, we aimed to identify novel genomic regions for the net CO2 assimilation rate (A) by combining genome-wide association study (GWAS) and the newly developed rapid closed gas exchange system MIC-100. Using three MIC-100 systems in the field at the vegetative stage, we measured A of 168 temperate japonica rice varieties with six replicates for three years. We found that the modern varieties exhibited higher A than the landraces, while there was no significant relationship between the release year and A among the modern varieties. Our GWAS scan revealed two major peaks located on chromosomes 4 and 8, which were repeatedly detected in the different experiments and in the generalized linear modelling approach. We suggest that high-throughput gas exchange measurements combined with GWAS is a reliable approach for understanding the genetic mechanisms underlying photosynthetic diversities in crop species.

Integrated genome-wide differentiation and association analyses identify causal genes underlying breeding-selected grain quality traits in japonica rice.

Improving grain quality is a primary objective in contemporary rice breeding. Japanese modern rice breeding has developed two different types of rice, eating and sake-brewing rice, with different grain characteristics, indicating the selection of variant gene alleles during the breeding process. Given the critical importance of promptly and efficiently identifying genes selected in past breeding for future molecular breeding, we conducted genome scans for divergence, genome-wide association studies, and map-based cloning. Consequently, we successfully identified two genes, OsMnS and OsWOX9D, both contributing to rice grain traits. OsMnS encodes a mannan synthase that increases the white core frequency in the endosperm, a desirable trait for sake brewing but decreases the grain appearance quality. OsWOX9D encodes a grass-specific homeobox-containing transcription factor, which enhances grain width for better sake brewing. Furthermore, haplotype analysis revealed that their defective alleles were selected in East Asia, but not Europe, during modern improvement. In addition, our analyses indicate that a reduction in grain mannan content during African rice domestication may also be caused a defective OsMnS allele due to breeding selection. This study not only reveals the delicate balance between grain appearance quality and nutrition in rice but also provides a new strategy for isolating causal genes underlying complex traits, based on the concept of "breeding-assisted genomics" in plants.

Independent control of organ number and distribution pattern in rice panicle.

As the determinants of yield products, rice panicle traits are important targets for breeding. Despite their importance in grain filling and subsequent yield productivity, knowledge on the organ distribution pattern in rice panicles is limited owing to the lack of objective evaluation methods. In this study, we developed a method for quantifying rice panicle organ distribution patterns. To validate our method for practical application in biology, we integrated this method into a quantitative trait locus (QTL) analysis and identified QTLs for panicle organ distribution patterns in rice. Interestingly, Grain number 1 (Gn1), a major QTL of organ number, was not identified as a QTL for distribution pattern, indicating that the number and distribution of panicle organs are independently controlled. This study provides insight into rice panicle organ distribution patterns that will help improve breeding targeting rice panicle architecture.

Effective use of legacy data in a genome-wide association studies improves the credibility of quantitative trait loci detection in rice.

Genome-wide association studies (GWASs) are used to detect quantitative trait loci (QTL) using genomic and phenotypic data as inputs. While genomic data are obtained with high throughput and low cost, obtaining phenotypic data requires a large amount of effort and time. In past breeding programs, researchers and breeders have conducted a large number of phenotypic surveys and accumulated results as legacy data. In this study, we conducted a GWAS using phenotypic data of temperate japonica rice (Oryza sativa) varieties from a public database. The GWAS using the legacy data detected several known agriculturally important genes, indicating reliability of the legacy data for GWAS. By comparing the GWAS using legacy data (L-GWAS) and a GWAS using phenotypic data that we measured (M-GWAS), we detected reliable QTL for agronomically important traits. These results suggest that an L-GWAS is a strong alternative to replicate tests to confirm the reproducibility of QTL detected by an M-GWAS. In addition, because legacy data have often been accumulated for many traits, it is possible to evaluate the pleiotropic effect of the QTL identified for the specific trait that we focused on with respect to various other traits. This study demonstrates the effectiveness of using legacy data for GWASs and proposes the use of legacy data to accelerate genomic breeding.

Genome-wide association study identifies a gene responsible for temperature-dependent rice germination.

Environment is an important determinant of agricultural productivity; therefore, crops have been bred with traits adapted to their environment. It is assumed that the physiology of seed germination is optimised for various climatic conditions. Here, to understand the genetic basis underlying seed germination, we conduct a genome-wide association study considering genotype-by-environment interactions on the germination rate of Japanese rice cultivars under different temperature conditions. We find that a 4 bp InDel in one of the 14-3-3 family genes, GF14h, preferentially changes the germination rate of rice under optimum temperature conditions. The GF14h protein constitutes a transcriptional regulatory module with a bZIP-type transcription factor, OREB1, and a florigen-like protein, MOTHER OF FT AND TFL 2, to control the germination rate by regulating abscisic acid (ABA)-responsive genes. The GF14h loss-of-function allele enhances ABA signalling and reduces the germination rate. This allele is found in rice varieties grown in the northern area and in modern cultivars of Japan and China, suggesting that it contributes to the geographical adaptation of rice. This study demonstrates the complicated molecular system involved in the regulation of seed germination in response to temperature, which has allowed rice to be grown in various geographical locations.

Evolutionary alterations in gene expression and enzymatic activities of gibberellin 3-oxidase 1 in Oryza.

Proper anther and pollen development are important for plant reproduction. The plant hormone gibberellin is important for anther development in rice, but its gametophytic functions remain largely unknown. Here, we report the functional and evolutionary analyses of rice gibberellin 3-oxidase 1 (OsGA3ox1), a gibberellin synthetic enzyme specifically expressed in the late developmental stages of anthers. Enzymatic and X-ray crystallography analyses reveal that OsGA3ox1 has a higher GA7 synthesis ratio than OsGA3ox2. In addition, we generate an osga3ox1 knockout mutant by genome editing and demonstrate the bioactive gibberellic acid synthesis by the OsGA3ox1 action during starch accumulation in pollen via invertase regulation. Furthermore, we analyze the evolution of Oryza GA3ox1s and reveal that their enzyme activity and gene expression have evolved in a way that is characteristic of the Oryza genus and contribute to their male reproduction ability.

Effects of irrigation on root growth and development of soybean: A 3-year sandy field experiment.

Increasing the water use efficiency of crops is an important agricultural goal closely related to the root system -the primary plant organ for water and nutrient acquisition. In an attempt to evaluate the response of root growth and development of soybean to water supply levels, 200 genotypes were grown in a sandy field for 3 years under irrigated and non-irrigated conditions, and 14 root traits together with shoot fresh weight and plant height were investigated. Three-way ANOVA revealed a significant effect of treatments and years on growth of plants, accounting for more than 80% of the total variability. The response of roots to irrigation was consistent over the years as most root traits were improved by irrigation. However, the actual values varied between years because the growth of plants was largely affected by the field microclimatic conditions (i.e., temperature, sunshine duration, and precipitation). Therefore, the best linear unbiased prediction values for each trait were calculated using the original data. Principal component analysis showed that most traits contributed to principal component (PC) 1, whereas average diameter, the ratio of thin and medium thickness root length to total root length contributed to PC2. Subsequently, we focused on selecting genotypes that exhibited significant improvements in root traits under irrigation than under non-irrigated conditions using the increment (I-index) and relative increment (RI-index) indices calculated for all traits. Finally, we screened for genotypes with high stability and root growth over the 3 years using the multi-trait selection index (MTSI).Six genotypes namely, GmJMC130, GmWMC178, GmJMC092, GmJMC068, GmWMC075, and GmJMC081 from the top 10% of genotypes scoring MTSI less than the selection threshold of 7.04 and 4.11 under irrigated and non-irrigated conditions, respectively, were selected. The selected genotypes have great potential for breeding cultivars with improved water usage abilities, meeting the goal of water-saving agriculture.

Potential of rice landraces with strong culms as genetic resources for improving lodging resistance against super typhoons.

It is generally believed that rice landraces with long culms are susceptible to lodging, and have not been utilized for breeding to improve lodging resistance. However, little is known about the structural culm strength of landraces and their beneficial genetic loci. Therefore, in this study, genome-wide association studies (GWAS) were performed using a rice population panel including Japanese rice landraces to identify beneficial loci associated with strong culms. As a result, the landraces were found to have higher structural culm strength and greater diversity than the breeding varieties. Genetic loci associated with strong culms were identified, and it was demonstrated that haplotypes with positive effects of those loci were present in a high proportion of these landraces. These results indicated that the utilization of the strong culm-associated loci present in Japanese rice landraces may further improve the lodging resistance of modern breeding varieties that have relied on semi-dwarfism.

Making the 'Green Revolution' Truly Green: Improving Crop Nitrogen Use Efficiency.

Traditional breeding for high-yielding crops has mainly relied on the widespread cultivation of gibberellin (GA)-deficient semi-dwarf varieties, as dwarfism increases lodging resistance and allows for high nitrogen use, resulting in high grain yield. Although the adoption of semi-dwarf varieties in rice and wheat breeding brought big success to the 'Green Revolution' in the 20th century, it consequently increased the demand for nitrogen-based fertilizer, which causes severe threat to ecosystems and sustainable agriculture. To make the 'Green Revolution' truly green, it is necessary to develop new varieties with high nitrogen use efficiency (NUE). Under this demand, research on NUE, mainly for rice, has made great strides in the last decade. This mini-review focuses on three aspects of recent epoch-making findings on rice breeding for high NUE. The first one on 'NUE genes related to GA signaling' shows how promising it is to improve NUE in semi-dwarf Green Revolution varieties. The second aspect centers around the nitrate transporter1.1B, NRT1.1B; studies have revealed a nutrient signaling pathway through the discovery of the nitrate-NRT1.1B-SPX4-NLP3 cascade. The last one is based on the recent finding that the teosinte branched1, cycloidea, proliferating cell factor (TCP)-domain protein 19 underlies the genomic basis of geographical adaptation to soil nitrogen; OsTCP19 regulates the expression of a key transacting factor, DLT/SMOS2, which participates in the signaling of four different phytohormones, GA, auxin, brassinosteroid and strigolactone. Collectively, these breakthrough findings represent a significant step toward breeding high-NUE rice in the future.

Whole-genome sequence diversity and association analysis of 198 soybean accessions in mini-core collections.

We performed whole-genome Illumina resequencing of 198 accessions to examine the genetic diversity and facilitate the use of soybean genetic resources and identified 10 million single nucleotide polymorphisms and 2.8 million small indels. Furthermore, PacBio resequencing of 10 accessions was performed, and a total of 2,033 structure variants were identified. Genetic diversity and structure analysis congregated the 198 accessions into three subgroups (Primitive, World, and Japan) and showed the possibility of a long and relatively isolated history of cultivated soybean in Japan. Additionally, the skewed regional distribution of variants in the genome, such as higher structural variations on the R gene clusters in the Japan group, suggested the possibility of selective sweeps during domestication or breeding. A genome-wide association study identified both known and novel causal variants on the genes controlling the flowering period. Novel candidate causal variants were also found on genes related to the seed coat colour by aligning together with Illumina and PacBio reads. The genomic sequences and variants obtained in this study have immense potential to provide information for soybean breeding and genetic studies that may uncover novel alleles or genes involved in agronomically important traits.

Evolution of GA Metabolic Enzymes in Land Plants.

Gibberellins (GAs) play key roles in various developmental processes in land plants. We studied the evolutionary trends of GA metabolic enzymes through a comprehensive homology search and phylogenetic analyses from bryophytes to angiosperms. Our analyses suggest that, in the process of evolution, plants were able to acquire GA metabolic enzymes in a stepwise manner and that the enzymes had rapidly diversified in angiosperms. As a good example of their rapid diversification, we focused on the GA-deactivating enzyme, GA 2-oxidase (GA2ox). Although the establishment of a GA system first occurred in lycophytes, its inactivation system mediated by GA2oxs was established at a much later time: the rise of gymnosperms and the rise of angiosperms through C19-GA2ox and C20-GA2ox development, respectively, as supported by the results of our direct examination of their enzymatic activities in vitro. Based on these comprehensive studies of GA metabolic enzymes, we discuss here that angiosperms rapidly developed a sophisticated system to delicately control the level of active GAs by increasing their copy numbers for their survival under different challenging environments.

Diverse panicle architecture results from various combinations of Prl5/GA20ox4 and Pbl6/APO1 alleles.

Panicle architecture directly affects crop productivity and is a key target of high-yield rice breeding. Panicle length strongly affects panicle architecture, but the underlying regulatory mechanisms are largely unknown. Here, we show that two quantitative trait loci (QTLs), PANICLE RACHIS LENGTH5 (Prl5) and PRIMARY BRANCH LENGTH6 (Pbl6), independently regulate panicle length in rice. Prl5 encodes a gibberellin biosynthesis enzyme, OsGA20ox4. The expression of Prl5 was higher in young panicles resulting in panicle rachis elongation. Pbl6 is identical to ABERRANT PANICLE ORGANIZATION 1 (APO1), encoding an F-box-containing protein. We found a novel function that higher expression of Pbl6 is responsible for primary branch elongation. RNA-seq analysis revealed that these two genes independently regulate panicle length at the level of gene expression. QTL pyramiding of both genes increased panicle length and productivity. By combining these two genes in various combinations, we designed numerous panicle architecture without trade-off relationship.

A common allosteric mechanism regulates homeostatic inactivation of auxin and gibberellin.

Allosteric regulation is protein activation by effector binding at a site other than the active site. Here, we show via X-ray structural analysis of gibberellin 2-oxidase 3 (GA2ox3), and auxin dioxygenase (DAO), that such a mechanism maintains hormonal homeostasis in plants. Both enzymes form multimers by interacting via GA4 and indole-3-acetic acid (IAA) at their binding interface. Via further functional analyses we reveal that multimerization of these enzymes gradually proceeds with increasing GA4 and IAA concentrations; multimerized enzymes have higher specific activities than monomer forms, a system that should favour the maintenance of homeostasis for these phytohormones. Molecular dynamic analysis suggests a possible mechanism underlying increased GA2ox3 activity by multimerization-GA4 in the interface of oligomerized GA2ox3s may be able to enter the active site with a low energy barrier. In summary, homeostatic systems for maintaining GA and IAA levels, based on a common allosteric mechanism, appear to have developed independently.

Genome-wide expression quantitative trait locus studies facilitate isolation of causal genes controlling panicle structure.

The morphology of rice (Oryza sativa L.) panicles is an important determinant of grain yield, and elucidation of the genetic control of panicle structure is very important for fulfilling the demand for high yield in breeding programs. In a quantitative trait locus (QTL) study using 82 backcross inbred lines (BILs) derived from Koshihikari and Habataki, 68 QTLs for 25 panicle morphological traits were identified. Gene expression profiling from inflorescence meristems of BILs was obtained. A combination of phenotypic QTL (pQTL) and expression QTL (eQTL) analysis revealed co-localization between pQTLs and eQTLs, consistent with significant correlations between phenotypic traits and gene expression levels. By combining pQTL and eQTL data, two genes were identified as controlling panicle structure: OsMADS18 modulates the average length of the primary rachis and OsFTL1 has pleiotropic effects on the total number of secondary rachides, number of grains per panicle, plant height and the length of flag leaves. Phenotypes were confirmed in RNA interference knocked-down plants and overexpressor lines. The combination of pQTL and eQTL analysis could facilitate identification of genes involved in rice panicle formation.